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SUPERSONIC UNSTALLED FLUTTER 

by J. J. Adamczyk, M. E. Goldstein, and M. J. Hartmann 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 

SUMMARY 

Recently two flutter analyses have been developed at NASA Lewis Research 
Center to predict the onset of supersonic ujistalled flutter of a cascade of two- 
dimensional airfoils. The first of these analyzes the onset of supersonic flutter 
at low levels of aerodynamic loading (i.e., backpressure), while the second 
examines the occurrence of supersonic flutter at moderate levels of aerodynamic 
loading. Both of these analyses are based on the linearized unsteady inviscid 
equations ol gas dynamics to model the flow field surrounding the cascade. The 
details of the development of the solution to each of these models have been pub- 
lished. The objective of the present papei' is to utilize these analyses in a para- 
metric study to show the effeets of cascade geometry, inlet Mach number, and 
backpressure on the onset of single and multi degree of freedom unstalled super- 
sonic flutter. Several of the results from this study are correlated against experi- 
mental qualitative observation to validate the models. 

INTRODUCTION 

The problem ol flutter has long plagued the development of high speed com- 
pressor fan stages. The solution to this problem is often costly both in terms of 
time and money. For this reason engine manufacturers as well as government 
agencies are currently supporting numerous research programs in an attempt to 
better define regions of flutter instability and solutions to the problem. To date 
their research activities have uncovered two regions of the operating map of fan 
stages where flutter can be encountered at high sjK'eds. These regions arc shown 
schematically on the performance map of a typical high speed fan stage in figure 1. 
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Region 1, the zone of moderate to high backpressure supersonic flutter can extend 
from the stall line of the fan down to its operating line. Experimental evidence of 
this flutter mode is presented in reference 1. 

Region 11, the zone of low backpressure supersonic flutter can extend from 
wide open discharge to slightly below the operating line of the stage. Numerous 
analytical as well as e.xperunental papers have appeared in the open literature 
which document the existence of this flutter mode (see refs. 2 to 4). Recently 
NASA Lewis Research Center has developed analytical analyses for both of these 
flutter zones. In reference 5 an mialyses is developed for the region of moderate 
to high backpressure flutter, while reference d deals with the problem of low 
backpressure supersonic flutter. Roth of these analyses are based on the linear- 
ized Luisteady inviscid two-dunensional equations of gas ilynamics to describe tne 
flow field surroujidbig an infinite cascade of oscillating tliin airfoils in a supersonic 
stream. The details of the mathematical development of these models will not be 
presented in the present paper. The objective of the present paper is to utilize 
these analyses in a parametric study to show the effects of cascade geometry, inlet 
Mach muiiber, reduced frequency, and backpressure on super.sonic flutter. Several 
of the results from this study are correlated against I'.xperimental qualitative obser- 
vation to \-alidate the models. 

MODEL 1'X)RMUIATK)N 

The present analyses represents an incremental annulus of a fan stage as an 
infinite two-dimensional cascade of thui airfoils. In both analj ses the steady 
relative flow approaching the cascade is assumed to be supersonic, with a subsonic 
axial velocity component, and satisfies the Kantrowitz unique incidence relationship. 
This flow configuration causes the weak oblique leading edge shock wave to propa- 
gate upstream of the cascade. At moderate to high pressures it is assumed that the 
steady pressure rise across the cascade is produced by a system of normal shock 
waves lying within the cascade passages. Downstream of the shock waves the 
steady flow is unilorm and subsonic. Figure 2 shows a sketch of this steady flow 
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MU'chaiucal damping any small unsup^HJitoil motion impaitod to tho cascade will 
gi*ow in time, causing the systetn to fail. 

For a single degree of freedom pitching or plunging motion the nondimensional 
work done by the gas stream on an airfoil in the cascaile is (see ref. 7) 

W = (V^,;rlm{C,^j} (1) 

for pitching motion aiul 

W = h^,;rlm[Cj (2) 

for plunging motion, where the symbol lm{ [ tlenotes the imaginary ^xirt of the 
bracketed quantity. The \-ariables appi*aring in these two equations are 

ami whieli are tlie eomplex moment eoeffieient about the pitching axes, 
the complex lift coefficient, the am|ditude of the pitching oscillation, and tlie 
amplitude of the plunging oscillation, respc'ct ively. There will be a tendeni'v 
for single cU'gret* of freedv)m flutter to occur in eitiu'r the* pitc'liing or plunging 
mode whenever the' imaginary [>art of the eomplex monuMit or lift eoeffieient 
lu'comes positive. I’or two- degree of frec'dom eoupU'd flutter the' nondimensional 
woi’k done by the g-as sti'eam on an airfoil in tlu' cascade is 

W - sin y Keal[ t h^^n- lm[ f'p + ^ y lni{ (2) 

where the symbol Keal{ J denotes the real part of the complex bracketed quantity 
and y is the phase angle at an instant in time bi'twecm the pitciiing and plunging 
motion of an airfoil. I iilike single degree of fluttei- whose onsi't is indepi'iident of 
the amplitude ol the aii’foil motion the onset of coupled flutter as shown by equa- 
tion (;5) is strongly dependent on the* ratio of the pitch to plunge amplitudes and their 
associated phase shiit y. 

Computations were performed based on the theory of refiM'enees f) and ti to 
determine the effect of backpressure, inlet Mach number, reduced frei}uency, 
cascade solidity, and stagger angle on tlu* nondimensional aerodynamic woi'k 
for single degree of freedom pitching and plunging motion and for coupled pitching 
and plunging motion. These results are shown plotted as a function of the phase 
angle between two adjacent airfoils (i.e., intc'rblade phase angle) at all times. 

II III III I • I ? I I I ) I > I ' 
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SINGLE DEGREE OF FREEDOM FLFTTER 

The geometry of the cascade assumed in this study is, solidity 1.3, stagger 
angle 60° mdess othenvise noted. Pitching axis location is assumed to lie at mid- 
chord. Figures 3 to 7 show graphs of the work per cycle for a cascade undergoing 
simple harmonic pitching motion at moderate to high backpressures. The in pas- 
sage normal shock wave is assumed to lie slightly upstream of the passage en- 
trance. Results are shown in figures 3 to 5 for inlet Mach niunbers of 1.2, 1.4, 
and 1.6, respectively. The parameter varied in each of these figures is the 
reduced frequency based on semichord from a value of 0.25 to 1.0. Recalling 
that positive work per cycle implies instability provided the effect of mechanical 
damping is neglected, figures 3 to 5 show' that for Mach numbers between 1.2 and 
1.6 single degree of freedom pitching (i.e. , torsional) flutter can exist for reduced 
frequencies less than 0.25 at moderate to high backpressures. At Mach numbers 
above 1.2 this mode of flutter exists for reduced frequencies above 0.25 but ceases 
to exist at reduced frequencies greater than 0.5. The effects of cascade solidity 
and stagger angle on torsional flutter at moderate to high backpressures are show'n 
in figures 6 and 7. These results were computed based on solidities of 1.2 and 1.4 
and stagger angles of 50° and 70°. The reduced frequency and inlet Mach number 
w'ere held constant at 0. 5 and 1. 4. Although figures 6 and 7 show no region of 
instability they do show that reducing the solidity and decreasing the stagger angle 
have a slight stabilizing effect on single degree of freedom torsional flutter at 
moderate to high backpressures. 

The nondimensional work per cycle for a cascade pitching about mid-chord at 
low backpressures is shown in figures » to 12. The cascade geometry assumed 
in the computation of these results is identical to that assumed for the previous 
results unless othenvise noted. Figures 8 to 10 show the effect of Mach number 
and reduced frequency. From these results it appears that single degree of freedom 
torsional flutter will not e.xist at low' backpressures at reduced frequencies in excess 
of 0.50. This reduced frequency limit is greater than the corresponding limit 
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established li\)m the pi*evious set of I’esults tor moderate to high backpressure 
cascade operation. This implies that a cascade operating with a finite backpressure 
across it will be less susceptible to torsional flutter than a cascade operating at 
low bacK pressures. Tliis baekpressuri/.ing phenomenon lias been observed in fans 
where it causes the torsional flutter boundary to bend back as the pressure ratio 
across the fan is raised (see fig. 1). I'lie etTi*ets of cascade solidity and stagger 
;ingle on the work per cycle for torsional oscillation at low backpressures are 
shown in figures 11 ;uid 12. I'he values of tlu> parameters used in these computa- 
tions are identical to those used previously in c-omputing the results in figiires (5 
and 7. These results like their I’ounterparts for moderate Ui high backpressure 
show that reducing the stagger angle and reducing the cascade solidity have a 
slight stabili/.ing effect on torsional flutter at low backpressures. 

The work per cycle for a cascade undergtiing simple harmonic plunging motion 
at moderate to high backpressures is shown in figures 12 to 17. Figures 12 to 15 
show' the effects of reduced fretiuency ami Mach number. These results imply that 
the work per cycle will remain negative over all values of interblade phase angle 
provided the reduced frequency is slightly greater than 0.20 for inlet Mach numbers 
up to 1.6. As the inlet Mach number is reduced from 1.6 the transition reduced 
frequency decreases to appro.vimately 0. 15 at Mach number 1.2. The effects of 
solidity and stagger angle on the work per cycle are shown in figures 16 luid 17. 
increasing the solidity slightly enhances the .stability of the system, while increas- 
ing the stagger has a destabilii'.ing effect. 

The results for plunging motion at low bacK pressure are not presented because 
they showed that the work per cycle always remained lu-gative Thus the theory of 
reference 6 predicts that flutter in a pu.'e plur.ging mode cannot occur at low back- 
pressures. However, the results presented in figures 12 to 17 which were based 
on the theory of reference 5 show that single degree of freedom bending fluttei- can 
occur at moderate to high backpressures. Thus it can be concluded tliat back- 
pressuring tends to induce bending flutter. K.\perimental evidence to support this 
observation is provided in referenct* 1. 
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COUPLED FLUTTER 

Caita (ref. 8) showed the existence of coupled flutter in compressor rotors. 
Unlike coupled flutter in fixed wing aircraft where the coupling between the bend- 
ing and torsional mode is due to the aerodynamic forces, the coupling in com- 
pressor rotors is caused by mechanical restraining forces associated with part 
span shrouds and flexible disks. Carta in his analysis assumed the motion of 
the rotor disk to be 


( . 27m r 
cut 

N p 


S = Ae 

where A is the amplitude of the motion, t is time, cu is the circular frequency, 
n is an integer, N is the number of rotor blades, p is the pitch chord ratio, and 
r the peripherial distance around the wheel. If the rotor blades are rigidly fixed 
to the deforming disk the plunging and pitching amplitude of a blade section at a 
given radial location are 
h^ = A cos 6 

2 7mA i tt/ 2 
Q! = e 

° Np 

respectively, where 6 is the local stagger angle of the blade sections. The non- 
dimensional aerodynamic work per cycle associated with this motion is 


W = Att cos Cj^ - 


27m 


Real C 


M 


pN cos 6 

where the mteger n can take on both positive and negative values corresponding to 
either a backward or forward traveling wave along the disk rim. The lift and mo- 
ment coefficients appearing in this equation arc to be evaluated at an interblade phase 
angle of 

27m 


a = 


N 


Computations were performed to assess the influence of the twist bend coupling 


ratio 
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€ _ 27m 

Np cos 6 

inlet Mach numbei* and reduced Ircquency on the aerodynamic work per cycle at 
low and high backpressures. The results lor moderate to high backpressures 
are shown plotted as a fujiction ol the harmonic number n for specified values 
of reduced frequency and inlet Mach number in figures 18 to 2o. The cascade 
solidity and stagger angle are 1.3 and GO^, respectively. These results show 
that for n negative and large the gas stream is supplying energy’ to the cascade, 
while for positive values of n the energy flow is generally from the cascade to 
the gas stream. Hence, the vibrational waves traveling in the direction of rota- 
tion appear to be less stable than those traveling in the opposite direction. For 
small absolute values of n and for x'cduced frequencies greater than 0.25 the 
work per cycle remains negative for .Mach numbers greater than 1.2. Since the 
twist bend coupling parameter c is directly proportional to n it is seen that an 
effective means of suppressing coupled flutter is by mechanically controlling the 
ratio of twist to bend in the lower order vibrational modes. These results also 
show that increasing reduced frequency lias a stabilizing effect on coupled flutter 
but not to the extent it does for single degree of freedom flutter. 

The effect of lowering the backpressure across the cascade on the work per 
cycle is shown in figures 21 to 23. The trends of these results are quite similar 
to those for moderate to high backpressure. Increasing reduced frequency and 
limiting the ratio of twist bend coupling suppresses coupled flutter at low back- 
pressure, as it did at high to moderate backpressures. From these results it 
appears that proper choice of these two parameters can prevent coupled flutter 
from occurring over the entire high speed operating angle of a fan or compressor 
stage . 

CONCLUSIONS 

The results of numerous calculation based on the theory of references 5 and 6 
has been presented to show the effect of inlet Mach number, cascade geometry, 

I [I1 lif til II llltlli 
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reduced frequency, and backpressure on the susceptibility of a cascade to single- 
and multi-degree of fxeedom flutter. It was s>^own that increasing reduced fre- 
quency and backpressure had a stabilizing effect on single degree of freedom 
torsional flutter. It was also shown that single degree of freedom bending flutter 
could occur at moderate to high backpressure if the reduced frequency was below 

0. 25. For a coupled vibrational mode it was shown that the occurrence of flutter 
could be prevented over the entire high speed operating range of a fan stage by 
mechanically controlling the ratio of twist bend coupling and reduced frequency. 
These stability trends are consistent with e.xperimental rig observations. 

A detailed experimental verification of the theories of references 5 and 6 
in necessary. One suitable approach would be to utilize a high speed wind tunnel 
with a linear cascade in which the airfoils are driven in a prescribed mode and 
operated over a range of back pressures. Measured surface pi'essure distribu- 
tions could then be correlated against predicted results to establish the model 
validity or the need for additional model refinements. 
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Figure 1. - Compressor performance and stability map. 
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Figure 3, - Work per cycle lor pitching motion atiout mid- 
chord dt moderate backpressure (Mach no. 1.2). 
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Figure 11. - The etiect of soliditv on the work per cycle tor 
pitching motion at low backpressure. 
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Figure I?. • The eftect of stagger angle on the work per 
cycle for pitching notion at low backpressure 
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Figure 13. - Work per cycle lor plunging motion at moder- 
ate backpressure iMach no. 1. Zi, 
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Figure 14. - Work per cycle for plunging motion at moder- 
ate backpressure (Mach no. 1.4). 



Figure 15. - Work per cycle for plunging motion at moder- 
ate backpressure (Mach no. 1.6). 



Figure 16. - The effect of solidity on the work per cycle for 
plunging motion at moderate backpressure. 
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Figure 17. - The effect of stagger angle on the '••ork per 
cycle at moderate backpressure. 


REDUCED FREQUENCY, 
k 

O 0.25 

□ .50 

A 1.00 



-30 ■ - i ! i 1 1 * 

-12 -8-40 48 12 

HARMONIC NUMBER 


Figure 18. - Work per cycle for coupled motion at moder- 
ate backpressure (Mach no. 1. 2i. 
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Figure 19. - Work per cycle (or coupled flutter at moderate 
backpressure (Mach no. 1. 4i. 
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Figure 20. - Work per cycle (or coupled motion at moder- 
ate backpressure (Mach no. 1. 6). 
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Figure 22. - Work per cycle for coupled motion al low back- 
pressure (Mach no. 1.4). 
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Figure 23. - Work per cycle for coupled motion at low back- 
pressure (Mach no. 1.6). 


